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Abstract. We construct local minimizers to the Ginzburg-Landau functional of supercon- 
ductivity whose number of vortices N is prescribed and blows up as the parameter e. inverse 
of the Ginzburg-Landau parameter k, tends to 0. We treat the case of N as large as |log e\, 
and a wide range of intensity of external magnetic field. The vortices of our solutions arrange 
themselves with uniform density over a subregion of the domain bounded by a "free boundary" 
determined via an obstacle problem, and asymptotically tend to minimize the "Coulombian 
renormalized energy" W introduced in |14j . The method, inspired by j22j . consists in min- 
imizing the energy over a suitable subset of the functional space, and in showing that the 
minimum is achieved in the interior of the subset. It also relies heavily on refined asymptotic 
estimates for the Ginzburg-Landau energy obtained in [14] . 



1. Introduction 

Let n C be a bounded and smooth domain. We are interested in finding local minimizers 
of the 2D Ginzburg-Landau functional of superconductivity 

(LI) Ge{u,A):=l- [ \{V-zA)u\'+\VxA-K,f+ ^^~J''J ^\ 

Here, the complex-valued function m is a pseudo wave function, the "order parameter" in 
physics, that indicates the local state of the superconductor. The material is in the supercon- 
ducting state in the regions where |n| ~ 1, and in the normal state where \u\ ~ 0. The zeros of 
the order parameter carrying a non-zero degree are called vortices. The parameter hex > is 
the intensity of the applied magnetic field, A is the induced magnetic potential, h := V x A is 
the induced magnetic field, and finally e is the inverse of the "Ginzburg-Landau parameter." 
For more background, we refer to standard physics textbooks, such as [211 [23], or to [19] and 
references therein. 

We are interested in the asymptotic regime e — )■ 0. This corresponds to the characteristic 
lengthscale of the vortices going to 0. 

The vortices of a critical point of flLip are conveniently studied via the vorticity measure 

A) ■=V X j{u, A) + V X A, 

where j{u, A) := {iu, Vau) denotes the superconducting current associated to a configuration 
and Va = V — iA. We also write := ^(me,/!^), for {u^yA^) a pair depending on e. In the 
limit e 0, the vorticity measure is well approximated by a sum of Dirac masses at the 
vortex centers, for a precise statement see the Jacobian estimate flL7p below. 

In what follows we write, for quantities and B,. depending on e, A,. < when there is a 
constant C independent of e such that A^ < CBs- We also write <^ B^ when lim^-^o ^ = 0- 
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The asymptotic study of critical points of f ll.ip . as e — ?■ for external fields /igx ^ ^2, has 
been carried out in a series of papers by Sandier and Serfaty, starting from [201 [21]; continuing 
in [16|, [TTJ [18]; the monograph [19], and more recently [H]. It was shown in [16] that for 
hex = o (I log e\) there is a critical value A* of A = lime_!.o ji^^ below which minimizers have 
no vortices. It was further shown in [17] that if hex — A* |log e\ < log |log e\ then the number of 
vortices of a minimizer remains bounded independent of e. For log |log e\ <^ hex — A*|log e\ <C 
|log e\, minimizers have vortices with 1 ^ ^ hex as shown in [19], Chap. 9, while for 
A > A* possibly = +oo, minimizers exhibit a number proportional to hex of vortices with 
uniform density in a subregion determined by an obstacle problem (see [18], [19] Chap. 7, and 

The behavior of critical points of fll.ip that are not necessarily minimizers of the energy, 
was the object of study of [18]. In [18] the vorticity measures fi{us,A^) associated to critical 
points {u^,Ai,) of fll.ip . normalized by the number of vortices, are shown to converge weakly 
to a measure /i. The measure /i has the form /i := —Ah + h, where /i is a "stationary function" 
for a free boundary problem (see f ll.2p below), and satisfies in some weak sense the relation 
/iV/i = 0. This expresses that the average force acting on the vortices (Vh) vanishes where 
the vortices are, i.e. on the support of /i, thus allowing them to be at equilibrium. This result 
was later proved under less restrictive assumptions in [19], Chapter 13, and included cases not 
previously covered. 

It is then natural to ask the converse, that is, given a measure fi associated to a stationary 
function h as above, can one find {ui;,A^) critical points of f ll.ip such that fi{ui;,Ai;), suit- 
ably normalized, converges weakly to //? We will answer this question partially, by building 
solutions to the Ginzburg-Landau equations which are local (but not global) minimizers of 
the energy, and have a prescribed and divergent (as e — 0) number of vortices A^. The same 
thing has already been accomplished for a bounded number of vortices in [22], and also in 
[19], Chap. 11 for A^ slowly diverging as e — 0. 

Note that the question of inverse problems, i.e. whether given a critical point of a limiting 
energy, one can find corresponding critical points of the original problem, has been addressed 
in a variety of contexts, and for Ginzburg-Landau in particular : see [H H] for the situation of a 
finite number of vortices in 2D, and [Bliniill] for lines of vortices in 3D. Our approach, which 
is not local inversion but rather local minimization, is more in the spirit of [TUl [TT] for vortex 
lines in 3D. All previous results however, concerned only finite numbers of vortices. Another 
related context where such a question is addressed is that of the AUen-Cahn equation, whose 
limiting energy is the perimeter functional, for results there see [H [71 [T3] . 

The stable solutions found in [22] and [19] Chap. 11, are shown to exist for a wide range 
of hex (they are in fact globally minimizing for only one value of hex)'- they exist for fields of 
intensity up to ~ where < |, and also for fields almost as small as constant. They thus 
form branches of solutions, each corresponding to a number A^ of vortices, with A^ bounded 

-2ri 

with £, or even unbounded, but subject to the restrictions A^^ < rjhex, hex < Ne^^^, where 

77 > is a small but fixed parameter. These restrictions imply in particular (cf. [19]) that 

1 

in any case A^ ^ (|log e\)^ and if A^ — 00, then hex ^ ^ for any a > 0, thus leaving 

open the question of existence of stable solutions with (|log e\)^ < A^ < hex vortices. We 
address this question and prove for the first time that branches of local minimizers with A^ 
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vortices continue to exist up to = 0{hex), with restrictions on and h^x specified below, 
see (fT3D-(fT:B]). 

The analysis in [22| [T9] consisted, roughly, in minimizing the energy Ge "among configura- 
tions with vortices" and exploiting the quantization of through a careful expansion of 
the minimal energy of a configuration in terms of its vortices. This expansion came with an 
error o(A^^). However, the method followed there only tolerates an error of o (|log e\) (hence at 
least the condition A^^ < |log e\) which thus limited the range of applicability of the approach. 
What we do here is to follow the same general idea of local minimization but this time over 
a class which is defined so as to exploit the results and methods of [H] to get a more precise 
energy expansion with an error of only o{N). In this way we can cover the range of A^ even 
comparable to |log e\, and fields as large as a power of ^. 

1.1. Statement of main result. In order to state the main result of the article in more 
accurate terms, we first proceed to introducing some auxiliary function. 

An Obstacle Problem. Let A^ be a positive number (typically an integer representing the 
number of vortices). Following [H], we denote by /i^ tv the minimizer of 

(1.2) i / \Vh\^ + \h-Kx\\ 

over the set of functions h that are equal to hex on the boundary and such that 

-Ah + h\= 2'kN. 

As long as 2'kN < hex\^\ (an assumption we will need to make), this minimizer h^^N exists 
and it is the solution to an obstacle problem (see Lemma 5.1 and Appendix A in [S]). There 
exists a coincidence set cJe^tv C Q and a number < < 1 (for us the normalized vortex 
density) such that 

(1.3) fis^N ■= -Ahs^N + hs^N = hex me,Arl^^^ and hex rris^N Ic^e.Af I = '2nN, 

where 1 denotes the characteristic function of a set and | ■ | its area. If Q is convex, then Us,n 
is also convex, by a result in [3]. To avoid technical difficulties, we will assume in the sequel 
that Q is convex. It is also checked in [H] that is bounded universally from below by a 
positive constant mo and that m^^Tv ^ \^e,N\ is increasing. Without loss of generality we can 
thus assume that 

(1.4) limmeAT =: m > mo > 0, 

where A^ may also depend on e. 

As in [13] (cf. Lemma 5.2), the function /le^r will serve to "split" the energy in a 
convenient way. The vorticity of the solutions we construct will ressemble fie^N, which is a 
constant over the subregion u^^n, hence we will have a uniform density of vortices there, and 
a vanishing one outside. 
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Range of fields and vorticity studied. We consider fields and vortex numbers satisfying tlie 
following conditions 

(1.5) K iV < min I c ^'"^' , K log ^ 



27r S\/he 
for constants < c < 1 and i^' > independent of e, and 

(1.6) hex < for some Sq < ^. 

Note that < log — could equivalently be written < /"^llog el (changing K if neces- 
sary). Let us also point out that (11. 5p implies hex ^ 1 and < rriQ < m < 1 (see (11.31) ). and 
thus 1 — rris^N is always bounded away from for e small enough, and /?.ex-(l — ^e,N) ^ 1- We 
will make use of these facts repeatedly during the proofs. 
Our main result is: 

Theorem 1.1. Let C 6e a convex smooth hounded set. Assume that N = N{e) a family 
of positive integers and hex = hex{s) satisfy conditions (II. 5p - (11.61) . Then, for e small enough 
(depending on the constants above), there exists a local minimizer {u^, A^) of (11.11) such that 
the following holds: 

(i) There exists a measure of the form 27t diSai, with di E Z and ai G Vl, such that 
(1-7) ||^(^., A,) - ye\\(^cl-^(n)Y < ^" 

for some a > 0, and we have 

(1.8) UeiVl) = 27lN, 

i.e. the total vorticity of is N. 

(ii) The asymptotics 

(1.9) G,{u,,A,) = + nN\og^= + NWm - 2nNhex{l - m,,;v) + o(iV), 

eWhp 



''ex 



hold as e — )■ 0, where Wm and are constants explicited in (II. lip and (I1.12p respec- 
tively, and m is as in (ll.4p . 

(iii) For any 1 < p < 2, 

'1-10) ||/i(Me, A^) - IJ,,^N\\w-i,P(n) < Cp\^, 

where Cp depends only on p. 

(iv) For 1 < p < 2, denoting by the probability measure on Lf^^(M^, M?) which is the push- 
forward of the normalized uniform measure on Ui^^jy by the map x H- V hexje ( V hex + 
■)),Q we have that, up to subsequence, Ps converges weakly to a translation invariant 
probability measure P on LfQ^(M^, M^) such thatP-a.e. j G Am andP-a.e. j minimizes 
W over Am, where W and Am o-i"^ defined just below. 



""^Here je is extended by outside fl. 
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As announced, this theorem proves the existence of locally minimizing solutions of Ginzburg- 
Landau with vortices. Moreover, their vortices tend to follow the average distribution 
/ig^AT, and their superconducting currents, after blow-up at scale -y/Zw, tend to minimize the 
"Coulombian renormalized energy" W, introduced in [11], whose definition we recall below. 

1.2. The renormalized energy W. For each m > 0, consider the class of currents j such 
that div j = and V x j = J2pev ^^^p — m for a discrete set V, and #(7^ fl B^) / is bounded 
for i? > 1. We call this class Am- 

For any family of sets {U/{}ij>o in R^, xur will be a family of associated positive cut-off 
functions satisfying 

iVxUfll < C, Supp(xufl) C \Jr, Xunix) = 1 if d{x, {VrY) > 1, 

where the constant C is independent of R. We specialize this to the family of balls Br of 
radius R and center the origin. 

Definition 1.2. W is defined, for j G Am, by 

= limsup .^^-^''^^^^ 
where for any positive function x, we denote 



\B, 



X) = lim I ^ / X bf + 7rlog?7^x(p) 

jR2\Up6pB(p,r,) 

For properties of W, we refer to [H]. It is proven there in particular that min^g^^ W is 
achieved and finite. We also set 7 



7 := lim 



- / + ^ + —dr- vrlni? 

2./ ^ 2 r 



constant introduced in [2], 



where / is such that f{r)e^^ is the unique (cf. P), up to a phase shift, degree one radial 
solution of 

- Am + (1 - \u\^)u = in R^. 

We finally let 

(1.11) := — min + 7. 



In [M], it is conjectured that the minimum of W over Am is achieved for currents corre- 
sponding to a configuration of points in a perfect hexagonal lattice (of suitable volume). A 
small hint to this is provided there by the fact that W is uniquely minimized, among perfect 
lattices of fixed volume, by the hexagonal lattice. In superconductivity, the hexagonal lattice 
also corresponds to the famous "Abrikosov lattice" , predicted from (11.11) and observed in ex- 
periments. It was proven in [H] that global minimizers of the Ginzburg-Landau functional 
(II. ip have currents which minimize W , after a suitable blow inside the free-boundary region 
that contains the vorticity. If the conjecture about the minimum of W is true, then one 
expects the vortices of global minimizers to arrange themselves in such an Abrikosov hexag- 
onal lattice. From item (iv) of Theorem 11.11 the exact same thing is expected of the lo cally 
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minimizing solutions found here. The only qualitative difference with the global minimizers 
is in the total number of vortices N, which manifests itself in a different /i^ tv hence a different 
local density of vortices and a different location for the free-boundary that encloses them. 

1.3. The local minimization method. The local minimization first relies on the same 
energy splitting introduced in [13], which we now recall. For any A, set 

A^:=A- V^h.^N 

where hs^N is the solution of the obstacle problem (II. 2p above. Let us also define accordingly 



(1 - 

J,^n{u,A):=^ I \{V -iA)u\'' + \V X A- fi.^^f + ^—^^dx, 



- [ |(V - iA)uf + \V X A- fi, 
2 Jn 

with /ie^Tv as in (11. Sp . the configuration independent quantity 

(1.12) Gf := \\h,^N - hexWli^n) + 27riV(m,,;v - l)hecc, 

and finally the function C, '■= hex — h^^N (implicitly depending on e and A^, note that it is 
nonnegative in Q and vanishes on dQ). We denote 

(1.13) ^max = ^^^^ = f^exi^ - ^s.n)- 

Lemma 1.3. For any {u,A), Gi.{u,A) can he decomposed as: 

(1.14) a(M,A) = Gf + J,,^(m,A)- / i^l{u,A^) + ]- [ {\uf-l)\Vh,,Nf. 

Jn ^ Jn 

Proof. This splitting formula was written with slightly different notation in [13]. Decomposi- 
tion (I1.14P follows easily from the relations: 

|(V - tA)uf = |(V - iAi)uf + \uf I V^/i,,jvf - 2j{u, Ai) ■ V^h.^N, 

and 

\VxA-hexf = \V X Ai + Ah^N - Km + h,M - Kxf 



\V X Ai- fis,Nf + \he,N - Kxf + 2(V X Ai- ^ie,N){he,N 



□ 



The last term in the right-hand side of (ll.Mp will be proven later to be o(l) for the config- 
urations of interest to us. 

As explained above for the case of smaller values of iV, we still seek our local minimizers 
as configurations that minimize among configurations with vortices. The question is to 
define (following [22]) a suitable set corresponding to this heuristic idea. This is achieved here 
by minimizing Ge over 

(1.15) 

S-f^ := {{u,A) e H\n,C) X H\n,R'^) s.t. J,Ar(M,Ai) < ttA^ log + A^VT^ + /3A^; 

h^i> 2tiNKx{1 - m,,^) - PN 
Jn 
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where /3 is chosen to be a small enough positive constant. The idea of this class is that it 
allows for at most vortices when /3 is small, it leaves enough room for a particular test 
function to be admissible (the precision with which we can compute its energy is o{N) while 
the error allowed is > N), and it can easily be related to the Ginzburg-Landau energy Gs 
via fll.Mp . As in [22| [T9] . the task is then, exploiting the quantization of A^, to show that 
miuc-s Ge is achieved in the interior of j^, thus providing a local minimizer of Gs, and to 

show it has the other desired properties. The class 5*^^ differs from those in [221 US] • Indeed, 
following these works it would be tempting to try minimizing over the set of configurations 
[u, A) such that: 



J,,^(m,A)- { nN log ^= + NW, 



m 



However, if nIu, Ai) = vrA^log 1 — + NWm — PN we cannot conclude that u has strictly 

' Sy Ilex 

less than N vortices, which is crucial to assure that such a configuration on the boundary of 
the set would have an energy that exceeds that of a local minimizer. If one replaces PN by 
/3 log / a similar problem occurs; in this scenario one is able to prove that if iu, A) is in the 

Ey hex 

corresponding lower boundary then it has at most — 1 vortices, but the —(3 log —k — in the 

£~\/ hex 

value of J(u, Ai) makes the energy of the hypothetical configuration possibly drop below the 
one in Proposition 13. ![ rendering that approach doomed. To overcome this difficulty (and this 
is a novelty of the paper) we proceed by bounding the total vorticity from below indirectly by 
requiring the elements of S'f ^ to satisfy 



e/i(u, Ai) > 2'kK,{1 - ms,N)N - f3N, 

n 

this implies there are at least vortices as long as there are no vortices with negative degree 
outside the free boundary u^^n- This is proven to be the case in Proposition 14. 1[ 

Finally, the question of treating the case of A^ larger than the order |log e\ remains open, 
we believe it would require much more precise estimates on the energy than those provided 
by [Hj that we use here, and this seems out of reach for the moment. 



The plan of the paper is as follows. In Section 2 we introduce notation and recall the 
Jacobian estimate. In Section 3 we construct an explicit configuration with A^ vortices uni- 
formly distributed according to fie,N and a superconducting current almost minimizing W, 
we check the test configuration belongs to the admissible class, and by computing its energy 
Gg, we obtain an almost optimal upper bound. In Section 4 we analyze general properties 
of the elements of 5*^^, in particular characterize their total vorticity as well as give lower 
bounds for the different components of the splitting of G^- We conclude with the proof of Theo- 
rem [TTT] in Section 5, by showing minimizers do not lie on the boundary of the admissible class. 
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2. Preliminaries 

We will mostly work in blown-up coordinates at the characteristic lengthscale l/\/h^. Let 
(2.1) e' := e^/h^ and x := ^/h^x. 

The domains fig and (we^at)' correspond to Vt and uj^^n in the blown-up coordinates x'. We 
transform the rest of the variables in a natural fashion, that is 



(2.2) u\x') = u{x), ^'{x') = ^(x) and A'{x') = ^=A(x) 

V Zip 



''ex 



for any magnetic potential A. 

The blown up currents and measures are defined in terms of these quantities in a clear way, 
also denoted with a prime. 

We define, for U G fie, the sets: 

(2.3) U := {x ene/d{x,U) <1} 

(2.4) U ■={xe njd{x, dU) > 1} 
and 

(2.5) U ■.= {xe Vte/d{x, dU) > 2}. 

We now recall the basic concentration estimate relating the Jacobian and the total vorticity, 
a la Jerrard-Soner 1121. 



Theorem 2.1. (cf. [H], Theorem 5) Assume Geiue^A^) < i^')^^ , where s < 1. Then there 
exists a measure v'^ that depends only on (not on A^) of the form = '^'^'Ylii^i^ai-, for 
G Z and ai G fig , with 



(2-6) ll/i' - ^'ell(c0^i(Q,)). < {ey/hex)^Ge{ue,Ae), 

and there exists a constant Ci > such that for any measurable set E C 

f~ e 

(2.7) \u'^\{E) < Ci/^ " 



I log I 



where 

(2.8) e,(u,v4) := IVa^IV |V X A| 



'1 - In 



2(e 



/^2 



3. Local minimization class and energy upper bound 

As mentioned, we follow the idea of [22], and Chapter 11 of [19], and seek for local minimizers 
of Gs among configurations with "iV vortices" by minimizing Ge over the class j^, cf. f ll.lSp . 
with /3 to be determined later. From now on we write [us, A^) an element in the Argmin of 
Ge in 5*^ ^. The fact that such a (u^, Ae) actually exists follows from classical arguments (a 
proof can be adapted from [22] or Chapter 11 in [IS]). The task is then to show a minimizer 
cannot lie on the boundary. 
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We first show that the classes S*^^ are non-empty and provide a priori estimates on the 
energy of a minimizer in those classes. 

Proposition 3.1. The classes are non-empty. Furthermore, if {ue,A^) is a minimizer 
of Gs over S^j^, then, as e ^ , 

(3.1) A,) < Gf + vriVllog e'\ + NW^ - 27tNK,{1 - m,,^) + o{N). 

Proof. The proof relies on computing the energy of a carefully built test configuration (u^, A^). 
The construction is carried out in Theorem 6 in [T3] where the argument makes use of the fact 
that in their case corresponds to the number of vortices of a global minimizer. We show 
here that that the construction still works for satisfying (11.51) and also make explicit why 
the test configuration belongs to jy, an issue not present in the global minimization setting 
of [H]. 

In Theorem 6 of [14], a test configuration {u^, A^) is built in the following way. 
Step 1: 

A discrete set of points Ve, Vg C cu^^tv and \Ve\ = N is chosen. A current jl is built that 
satisfies 
(3.2) 



1 / 1 



lim sup — . r - 

The current jl is constructed under the assumption that Q is convex (Proposition 4.2 and 
Corollary 7.4 in [Hj). In addition it is necessary that the coincidence set Ue,N be well contained 
inside Q. More precisely one must guarantee 



(3.3) dist {Ue,N, fi"") > {me,Nhex) 2 

In [13] an explicit lower bound for dist {uJe,Ni ^'^) was proven in the case where is the 
number of vortices corresponding to global minimizers of G^. In our case we recall that (11.51) 
implies m < 1 and therefore -^1 — m^^jq is bounded from below by a positive number. From 
Appendix A in [T3], we know that 



(3.4) ||V/l£,Ar||oo < /iexA/l - ?Tie,Ar. 

Then since h^^N is equal to nie^N hex on uje,N and to hex on the boundary, ^ = hex — h^^N and 
(I1.13p . one has: 

dist {Ue.N, ^'') > u^^"'''\ — ^ a/1 - m^.AT. 
II Vfte,Ar||^oo 

On the other hand, m^^N > niQ > 0, 1 — nis^N is bounded below, and hex ^ 1, therefore 

/ ^ 1 

V '^e,N Ilex 

which proves (13. 3 p in our case. 
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Step 2: 

The functions and A,, are then built satisfying 
(3.5) luel = 1 on \ Upev,B{p, Me) and Up^v, B{p, Me) C Ue,N- 



(3.6) A, := Ai,, + V^/i,,Ar, 

where Ai^^ is chosen to satisfy V x Ai^^ = /i^^Tv and 



(3.7) - 
for all p EVe- 



B{p,Me) 



2 fi-ln; 



|2\2 



2^2 



7rlogM + 7 + OA/(l) + o,(l), 



(3.8) 



B(p,riime,N h^^)-^'^\B{p,Me) 



2 fl-IC 
+ 



|2\2 



2^2 



< TT log 



Me 



+ o{v). 



fl{Ue, Ai^e) = / fJ'e,N = '^nN, and /i(Ue, y4i,e) = OU (w^, 



(3.9) 

'n Jn 
From this it easily follows that (tTe, A^) G S^jv- Indeed 



1 



n 



1 

2./n 



V X - 



2 {1-\U, 



|2n2 



2^2 



2^ (1- Im; 



|2n2 



2e2 



Above, the second term disappears because of (13.61) . 

The identity |Vam|^ = |V \u\^ + \u'f implies, together with (13. 5p that 



= ISr onfi\Up6P,S(p,M,). 
Therefore, combining (13.21) . (13. 7p and (13. 8p . one gets after taking M — )■ oo, next ?7 — )■ 0, 
(3.10) Js^(n;,Ar^) < 7riV|loge| + 27riV min IV + tt log ^ +iV7 + o(iV). 

By assumption me^Ar — )■ m, therefore llogm^ Ar| = A^ |logm| + o{N). According to (1.17) and 
(1.18) of [H], by scaling, one has 

min W = — min W -\ — log m, 

^1 m Am 4 



Thus, by definition of Wm and e', (I3.10p is equal to 

(3.11) Je,N{ue,A^e) < 71 N \\oge' \ + NW^ + o{N) . 
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This shows that {us^^i,e) satisfies the upper bound in the definition of j^. To see it also 
satisfies the lower bound, we see by (13. 9p that 



(3.12) = ^rnax / ^l{Us, M,s) = ^max / fJ's,N = 27rNhex{l -rrie^N)- 

Jn Jn 

By fIXTTD and flXT^ we obtain that («;, A,) e S'f^^. 

To conclude, from the splitting formula (11.141) . we see that 

GMX) = + Je,N{^e,A^e) - I i^i{Ue.A^e) + \ I - 1) | V/l,,;v |' • 

The last term is identically zero because |V/ie.Ar| = on uJs^Ni while \ue\ = 1 on [us^nY by 
(13. 5p . Hence, by (13. lip and (I3.12p again, we see that 

Ge{ue, Ae) < Gf + vriV|log e'\ + NW^ - 2t:NK,{1 - m,,^) + o(iV), 
and the proposition is proved. □ 

In the next section, the lower bound component of Theorem 11.11 requires a very precise 
approximation of the vorticity measure by Ve- We also need a better control on the rightmost 
term in the splitting formula (I1.14p . 

Lemma 3.2. Let (m^, A^) he a minimizing pair on S^^. Assume N and hex satisfy (II. 5p and 
fll.6l) . Then there exists an a > such that 



(3-13) ll/i' - ^'e\\(^C°-\n,))* = 

and 

(3.14) / (1 - \Ue\^) |V/l,,^|' = 0(1). 

Proof. It is easy to see that < /i^^, while the rest of the positive terms in (13. ip are o{h1^) 
Therefore 

Ge{Ue,Ae)<hl^. 

From (12. 6 p and (II. 6p we can deduce 



(3.15) ll/i' - ^'e\\cr(n.)y < ehl ■ h^ < e^t^". 

The fact that this is o(e°) for some a > follows from (II. 6p . Thus, we have proved (I3.13p . 
and (I1.7P also follows. For the second assertion, Cauchy-Schwarz yields 



1 

J^{\Uef - 1) \Vh,,Nf < (^(1 - Kff^ ' ||V/l,,^||^ |1^|^ < y/e^G,{Ue,A,)\\Vhe, 



2 

A^lloo 
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Using (13 ■4p and the bound h^^ > Gs{u^, A^), the above leads to 



y"(l - 1^,1^) \Vh,M' < ehl < = 0(1), 



by ([ID. 

□ 

4. Further Properties of minimizers and lower bound for Gsiu£,Ae) 

The next proposition, which is the crucial part of the proof, shows that the elements of 
have vortices and also that the upper bound in (13. ip is sharp. 

At this point, it is convenient to introduce the notion of "vorticity adapted" family of 
functions. To that end, for a function ip^, let u^" := |x G S-t. ipsi^) = ||log ^'1} • We will 
restrict to considering families {ipe}e>o such that 

(4.1) \\VA\\oo = o{\\oge'\). 

Now, define the class: 

VA := I {V^£}e>o such that T/;^ : fie ]R,^£ satisfies (jlTT]), w'''" D (We_Ar)'; 

< ^pe{x) < ^|log e'\ if X e fie \ w'^- = on dQe 

Abusing notation, we sometimes call ip;; a vorticity adapted function when {ipe}e>o is a vorticity 
adapted family. This notion will help us in making use of a family of estimates available in 
[13] in a transparent way. The proof in [13] of such estimates, is presented for a particular 
vorticity adapted function, while here we use the estimates in their full generality. Ultimately, 
these estimates, together with the quantization of the vorticity measure, will prove that S*^^ 
consists only of functions with A^ essential vortices. 

In what follows, a + or — superscript denotes the positive resp. negative part of a function 
or measure. 



Proposition 4.1. Let N and hex be sequences satisfying (11.51) and (II. 6p . Assume {u'^,A'^) G 
S^j^. Then, fore small enough, we have: 

(i) There are no vortices with negative associated degree outside the free boundary, i.e.: 

(ii) The total vorticity is N, in other words: 

z/e(fi) = 271 N. 

Cm) 

^fi{u', Al) < 27lNhex{l - ms,N) + o(l). 



(ivj 

Je,N{u', Al) > TxN \\oge'\ + NW^ + o{N). 
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Proof. Recall that (m^, A'^) G means that 

(4.2) JeAu', Al) < ttN \\oge'\ + NWm + f3N, 
and 

(4.3) / e/i(w^ A\) > 27riV/ie..(l - ' I^N. 

Jn 

The proof relies on lower bounds on quantities of the form 

J,,n{u',AI)- / ^,{./hZ-)ij{u',Al), 
Jn 

where {tp£}e>o G VA. 

Step 1: First we claim we may assume that {u^, A"^) satisfies 

(4.4) div j{u\ A') =0. 

This assumption is justified by the following: we claim there is a Aq such that div j(u'^, A^) = 
and 

J.A^^AD- / ilj,{^/hr,-)fi{u',Al) > J,^n{u^AI,)- / tlje{Vhr.-)l^{u\Al^) + o{l), 
Jn Jn 

holds for any {■ipe}e>o £ VA Indeed, letting Aq be such that G'e(n^, Aq) = mmA^m(n,R'2) Ge{u^, A) 
(such a minimizer is easily seen to exist), by minimality {u^, Aq) satisfies the second Ginzburg- 
Landau equations, that is — V"'"(V x Aq) = (m^, V^g^i^), whence j{u'^,Al) is divergence free. 

By definition of Aq, Gs{u^, Aq) < Ge{u'^ , A'^) < /i^^. As in the proof of Lemma |3T2] this upper 
bound implies 

/ (1 - \U'\^) \Whe,N? = oil). 

Jn 

We can use the splitting formula f ll.l4p to conclude: 

Gf + J,,^(«^Af)- /^M^^D > Gf + 4^(n^A^J- /e/i(«^Ay + 0(1). 

Jn Jn 
From this, we can assert that 

47v(u^Af)- / ^/^(0^■)/^(M^A^) > J,,jv(m',^o,i) - / V^e(v^Xw',^o,i) 
Jn Jn 

[U^/hZ-) - i) Hn\A\) - ^i{u',Al,)) + 0(1). 
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If the field hex is not too large, the Jacobian estimate (12 ■6p implies that the last term is o(l). 
To see this, recall that the measure z/^ in (12.61) does not depend on A,, and therefore: 



v^(v^-)-e) (/i(w^Ag,l)-/i(w^AD) 



< 




■ max < 








oo 1 





< 



UlL + II V^lloo) ■ (v^Ge(«^ A^)) < {he. + d|log e\) hi = o(l) 



thanks to (II. 6p . This proves the claim. 

Step 2: Lower bound with vorticity adapted functions. We switch to blown-up coor- 



dinates in order to appeal to the proof of Proposition 6.1 in [13]. Recalling that the Jacobian 
estimate allows us to approximate fi' by z/^, we can replace bounding below (14. 4p by bounding 
below 



But it was proved in [H] that there exists a function g'^ depending only on {{u^)', (A^)') 
such that for any vorticity adapted function ip^, and for some constant Ci > 0, 



t/ J \ I O I / 1 

(4.5) +g,,{{Cj;:N)')+o{\{uJe,N)'\), 
where is as in (12. 8p . 

The argument can be found in (6.28)-(6.35) of [H] where C,'^ can be replaced by ipe'-, indeed 
it can be readily checked that the only fact about Ce that was used in the proof was that it 
was vorticity adapted. Note that here we also use the fact that div j{u^,A'^) = 0, obtained 
in Step 1. 

We also record the foUowings facts from [TU [T5] 

1. There exists a family of disjoint balls of bounded total radius such that the 

support of u'^ is contained in Ujg/i?j, a universal constant C2 and functions {f?^'} such 
that if u'^{B,) < 

(4.6) ^ g^,^ > I (|log e'\ - C2) W,{B,)\ > ^|log e'\ \u,{B,)\ 

(The family {Bi}i^j is built in [13], where (14. 6 p is proved in Proposition 3.2) 
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2. There exists a universal constant C3 such that 



(4.7) {9e'yix)>l\ Yl {9s'i^)-C3)lB. 



{i,u'^{B,)<0} 



(This is a consequence of the definition of Qs' in [14], combined with (3.22) in [15] and 

dH)). 

3. Let P be the probabihty measure in item (iv) of Theorem II .11 Because (m^, A^) satisfies 
the bounds in the definition of j^, and since the field h^x is not too large (cf. fll.6p ). 
the conclusions of Proposition 6.3 in [13] hold in our case with the density m instead 
of rrix- Assumption fl4.4p then allows us to assert (recalling the definitions (11. lip and 

) g'si^N) ^^{^j ^0') dPU) + 7 + 0(1)^ + o(|(u;.,7v)'|) > NW^ + o(iV), 
since, in view of (11.31) 



(4.9) \{u,,n)'\ = K. = ^<^ = 0{N). 

me,N mo 

Using (14.50 together with (14. 80 - 04. 9p . we have thus found that for any e VA, 

(4.10) J'sAi^n'. (Air) -J^^^du:>^J^ (l - + NW^ + o(iV). 



Step 3: All vortices of are contained in a neighborhood of (wg at)'. 

We will quite often need to use the fact that for any measurable set A, 
(4.11) u'^{A) e 27rZ, while \{u'^)+{A)\ , {A)\ , (A) e 2m. 

In what remains, a is set to be the constant 

I log e'\ 



a 



2hpJl — m 



e,N) 



Note that a depends on e and but we omit to indicate this dependence. Also, from here 
on, all the primed quantities denote the blown up versions of the unprimed ones, according to 
the transformations (12.11) - (12. 2p . We begin by noticing that a^' e VA. Indeed, since hex ^ 1, 
and m = lim£_5.o m^^N < 1 one can readily see, using (13.40 and the fact that hexi^ — ms^N) ^ 1, 
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that 
(4,12) 

\nex{^ - m^^N) J y/hex \ [hex) - m^.^N) 



oo 



<o( P^^^^ I Kx^/l-m,N<o( ^ '^"^^ ^ 1 =o(|log£'|). 

\{hex)^{l - me,N) J \Vhexil - m^,N) 

By (112]), (SJ]) and definition of a, we find that 
(4.13) JiJ{un',{Air)-a f eV((«7, (^i)') 



< 7riV|log e'\ + iVW„ + /SiV - a{2nNKxil - me,N) - l3N) 

= (1 + a)l3N + NWm. 

Conversely, the bound (14.101) apphed to ipe = together with (14.131) will allow us to rule 
out the presence of vortices with negative degrees where ^' is small. In turn, we shall see that 
this implies an upper bound on the total vorticity Uei^e)- 

As announced, we apply (14.101) to ip^ = a^' which is indeed vorticity adapted, combining 
with the Jacobian estimate (12. 6p . (14.131) and (14.101) we obtain that 



{l + a)PN + NW^ > J',^^{{u^)\{A{)')-a ^^ {{u^)\{A\y) 

> j;^((«7,(AD')-«^ e'c?^: + 0(liv«)llooll/i'-^:il(4>.i(^^))- 

(4.14) - IL + + + 

The last term is o(l) because Gs{u'^,A'^) < h^^ and hex satisfies (11.61) . Thus, we obtain 

For any Z > small, we introduce the following neighborhood of dfl: 



:= <^ x/e{x) < 1 ^ hexil - m,, 



TT 



where Ci and K are the constants in (12.71) and (11.51) respectively. From the definition of a, it 
follows that 



2<- ^ 6C,(i^+«g^)/3 _^ 

|loge'| ~ TT 
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Using this together with fl2.7p . we obtain (cf. (12 ■4p ) 



2 — ) 

(4.16) > ^ ^— (f//^) lloge'l . 

TT 

On the other hand, by (O]), < K\\oge'\ and jj^^ < 2n{iSi,M) ' combining f05|) 



les 



and (I4.16p . we get KM < 1 + o < 2 (for e small enough), which impli 

(4.17) H\ m = 0, 

in light of flCTl) . 

Step 3: There are no vortices with negative degree outside (we,Ar)' 

We turn to the proof of item (i). We invoke ( 14. 5 p applied to a^' and (I4.8p -( l4y9|) again to 
get 

j;^((n^)', (AIY) - [ ai'du'^ > ^(^,,)+((^T) + NW^m + o(iV), 
which, arguing as in (I4.14p yields 

(4.18) (1 + «)/3iV + o(iV) > ^{g^,)+{{cJZN'r)- 

Recall that from ([LS]) we have a(3N < (3 \\oge'\ (^^) < (3 \\oge'\ hence (1 + a)(3N + o{N) < 

/3|l0g£'|+0(|l0g£'|). 

On the other hand, from (14.61) and (14. 7p one can bound {gi;')~^{{uJe,N'y) from below by 
(letting C denote a generic positive constant) 

{ge'Vii^N'r) > c\\og e'\ \u;{{^^r)\ - cj-^^^^^^^f^ i^^i ^ ^'1 • 

Going back to (I4.18p . this implies /3+o(l) > C [(u^'Y) | , where C is a positive constant 
independent of (3 and e. Since G 27rN, we may choose /3 small so that for any e small 
enough, it imposes 

(4.19) {i^N'r)\ = 0- 
Since {uJ^'y contains {u'^^j^y, we have thus proved item (i). 

Step 5: The total vorticity of n"^ is 



We proceed, once more, applying (I4.10p only this time to a vorticity adapted function 
such that w'^^ ^ fie \ f/^. Such a function exists. Indeed, by (I4A2|) . for x G ^(fieV?/^) 
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we have < o{\\og e'\)d{x,dQe) while = ^(1 - 0{l)P)\\og e'\. It follows that 

(i(5(f2e\f/^), 3> 1 as £ — )■ (if /3 is a small enough constant), and thus (i((9(f2e\f/^), Sfi^) ^ 

1. One can therefore build a family of functions ipe which are equal to ||log e'\ in Q^\Uf^ = 

Qe\U^, which vanish on dQs and satisfy ||V^/'e||oo = o(|log e'\), i.e. belongs to VA. 
Applying fl4.10p to that family yields 

(A^i)') - [ Adiy',>NWm + o{N). 

We can then appeal to f l3.13p again to obtain 



> (ADO - / ^^du: + Oi\\V4^e\\^e^)>NW,n + o{N). 



Since {u^, A^) E 5'f ^, we must have 

(4.20) 7TN\\oge'\+m>Je,Niiu'y,iAiy)-NW^> [ ^edu', + o{N). 
Now, by (14.171) and definition of ips 

(4.21) / dv'^ = [ dv'^ =1^1 llog^'l dv'^ = \ lloge'l \ U^). 

Plugging this into fl4.20p leads us to 

B N 1 > ( N 



IT |log £'1 27r \ |log e' 

This implies that for e small enough, 

(4.22) N > ^(Mi^ = 

27r 2tt 

because of (14. lip and (I4.17p . We now show that there is actually equality in (I4.22p . To that 
end we first prove 

' ev((«7,(^iy)<^.(a)Cax+o(i)- 



This is not immediately clear since ^ is not constant and is not necessarily equal to z/+. 
However, from (11. 6p and (II. 7p . we have 



[eKiun',iAiy)= I i'dv', + 0{e'^hi:'\Vi'\J= li'dv'^+ ! i' dv'^ + o{e-) . 
On the other hand, from (14.191) . 
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Also, 

because ^' = ^max on (cJ^'). Collecting these estimates we are led to 

(4.23) / eV((«^)', {Al)') < W {ly'sii^N) + (z/:) + ((u^')')) + o{l) = U.ysi^e) + 0(1) 

where the last equality follows from (14.191) . Consequently i^'^i^e) > 27r(A^ — 1), for otherwise 
we would have (cf. (I1.13P ) 

/ ey((^')', (^i)') < 27r(iV - 1)^^,. + o(l) = 27r(iV - - m,,^) + o(l) 

and the right-hand side is strictly smaller than 2TTNhf.x{l — ttl^^n) — 13 N, if /3 is small enough, 
by virtue of (11. 5p . This contradicts {u'^,Al) G 5"^^. Therefore, ^'^(fig) > 2ttN (recall i^'^i^e) ^ 
27rZ.) Thus, with (I4.22p . we have proved item (ii). 



Step 6: Proof of (iii) and (iv) 

Item (iii) follows readily from (]4.23p since we already know that i^'^i^e) = ^.ttN. For item 
(iv), note that from (gTTjl . (Km and (lOTi) we deduce 

J;^((m^)', (AIY) > [ A du', + NW^ + o{N) 

= ^ \\oge'\ \ UP) + NWm + o{N) = irN \\oge'\ + NW^ + o{N). 
The proposition is proved. □ 



5. Conclusion 

In this section we finish the proof of our main result. 
Proof of Theoremll.il 

As mentioned earlier, a minimizer in S"^^ exists, so to show it actually corresponds to a 
local minimizer, it only remains to prove (m^, A^) ^ dS^ n- 
We decompose dS^j^ = (dS^^j^)'^ U(^'S'f at)^ where 

{dS^^^)+ := {(n, A)eH^x such that: Je,iv(«, ^i) = vrA^ |loge'| + NW^ + PN] 

and 

{dSl^Y := {{u, A)eH^ X such that: [ ^ii{u, A^) = 2TiNhe^{l - me,N) - PN). 
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Since (m^, A^) is a minimizer in S'^^y, we may apply Lemma 13^ and Proposition 14.11 By (11 .14^ 
and fl3.14p . we have 

(5.1) G,(m„ a,) = + J,,^(m„ - / e/iK, Ai,,) + 0(1). 

If {us,As) G {dS^ j^)^ then by (iii) in Proposition 14. H we have 

Geiu,, A,) > Gf + vriV |log£'| + NWm. + /3iV - 27rNh,,{l - m.^^) + o(iV). 
If instead, {ue,As) G {dS^j^)^, invoking (iv) in Proposition 14. II yields 

G,(m„ a,) > + ttN \\oge'\ + - (27riV/ie.(l - m,,^) - /3iV) + o(iV). 

In either case there is a contradiction with the upper bound (13. ip . so necessarily (u^, A^) ^ 
dS^ Pf, whence {u^, A^) is a local minimizer. 

From (15. ip . (iii) and (iv) in Proposition 14. H we must have 

(5.2) Geiue, A,) > Gf + vriV |log£'| + NW^ - 2tiNK,{1 - m,,;v) + o{N). 

Thus, the asymptotic expansion (II. 9p follows from (15. 2p and (13. ip . (II. 7p . (ll.Sp were proved in 
Lemma 13.21 and Proposition 14.11 respectively. 

We prove item (iii). It essentially follows the proof of (1.34) in [13]. First, arguing as in 
item 2 of Proposition 6.1 in [T3], replacing F'^ with J'^ ^ — a ^'fi', we can sandwich the term 

|V X Ai^s — /^e.ATp in Je{us,Ai^s) between the upper and lower bound of (I4.14p and hence 
control it by 0{N), i.e. 

(5.3) ||V X < TiV 

Using the same item of Proposition 6.1 in [T3] also allows to control, for 1 < p < 2, 

/ |j(<,Ayr<G,fj;^«,Ay-a / ^f^' + Kn\) <c{N+\uij). 

Rescaling we find that for any 1 < p < 2, 

(5.4) \\j{us,A^,Mmn) < ht' {N'/^ + H^M'^^) < hth'/^ < iV^ 

thanks to (14. 9p . then (II. 5p giving < hex- On the other hand, by direct calculation we have 

fl{Ue, Ae) - fl{Ue, Ai,^) = V X ((1 - \Ue\^)V^he,N) 

from (ODand ([M]) 

11(1 - \U,\') IVhML^ < V^e- \^K,nL ^^hl = o(l), 

from which one obtains 

(5.5) W^Ue, Ai^e) - Ku^, Ae)\\y^r-^,p = o(l). 

We may now write, using the triangle inequality 

||/i(ne. As) - fie,N\\w-^'P ^ ^'^) ^ /^(^e' II HZ-Lf + ll'"^:^ " V X Ai^e\\w-hP 

+ \\fi{Ue, Ai,,) - V X < N 
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where we have used (15.31) (together with an embedding inequahty), (I5.5p and the fact that 
li{us, g) — V X e = V X j(u£, Ai^^) combined with (15 ■4p . Consequently, item (iii) is proved. 
We turn to the proof of (iv). Since a^' G VA, (14. 5p . (14. 8 p and (ll.9p allow us to assert 

NWm + o{N) >N{^j WU) dP{j) + 7 + 0(1)^ 



which in turns translates into 



Wm>—[ W{j)dP{j)+J. 

m J 



To conclude, since P-a.e. j G Am (as given by [H]) this implies P-a.e. j minimizes W over 
Am- The result follows. 

□ 
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